Magnetic Frustration in a Mn Honeycomb Lattice Induced by Mn-O-O-Mn Pathways 



H. Wadati, x ' 2 'B K. Kate. 3 Y. Wakisaka, 3 T. Sudayama, 3 D. G. Hawthorn, 4 T. Z. Regier, 5 
N. Onishi, 6 M. Azuma, 6 | Y. Shimakawa, 6 T. Mizokawa, 3 A. Tanaka, 7 and G. A. Sawatzky 2 

1 Department of Applied Physics and Quantum-Phase Electronics Center (QPEC), 
University of Tokyo, Hongo, Tokyo 113-8656, Japan 
Department of Physics and Astronomy, University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada 
'' Department of Physics and Department of Complexity Science and Engineering, 
University of Tokyo, Kashiwa, Chiba 277-8561, Japan 
''Department of Physics and Astronomy, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada 
5 Canadian Light Source, University of Saskatchewan, Saskatoon, Saskatchewan S7N 0X4, Canada 
6 Institute for Chemical Research, Kyoto University, Uji, Kyoto 611-0011, Japan 
7 Department of Quantum Matters, ADSM, Hiroshima University, Hiroshima 739-8530, Japan 

(Dated: January 17, 2011) 

We investigated the electronic structure of layered Mn oxide Bi3Mn40i2(NOs) with a Mn hon- 
eycomb lattice by x-ray absorption spectroscopy. The valence of Mn was determined to be 4+ 
with a small charge-transfer energy. We estimated the values of superexchange interactions up to 
the fourth nearest neighbors (Ji, J2, J3, and J4) by unrestricted Hartree-Fock calculations and a 
perturbation method. We found that the absolute values of Ji through J4 are similar with pos- 
itive (antiferromagnetic) Ji and J4, and negative (ferromagnetic) J2 and J3, due to Mn-O-O-Mn 
pathways activated by the smallness of charge-transfer energy. The negative J3 provides magnetic 
frustration in the honeycomb lattice to prevent long-range ordering. 

PACS numbers: 71.30.+h, 71.28.+d, 79. 60. Dp, 73.61.-r 



Since the resonating valence bond state in geometri- 
cally frustrated magnets has been proposed by Anderson 
[l|, spin-disordered ground states in Mott insulators on 
frustrated lattices have been attracting great interest in 
condensed-matter physics. The exchange interaction J in 
a Mott insulator is roughly given by —2t 2 /E g , where t is 
the transfer integral between the two localized orbitals 
and Eg is the excitation energy across the Mott gap. 
In Mott insulators on frustrated lattices, spin-disordered 
systems including organic and inorganic materials 
all have relatively small E g , suggesting that the smallness 
of Eg or the closeness to the Mott transition would be 
important to realize the spin-disordered ground states. 

Various insulating transition-metal oxides are known 
as Mott insulators and can be classified into (i) the 
Mott-Hubbard type insulators where the Mott gap E g 
is mainly determined by the Coulomb interaction U be- 
tween the transition-metal d electrons and (ii) the charge- 
transfer type insulators where E g is determined by the 
charge-transfer energy A from the oxygen p state to the 
transition- metal d state 0]. Therefore, the smallness of 
E g can be obtained in transition-metal oxides with small 
U or small A. In the small U case, theoretical stud- 
ies on triangular-lattice Hubbard models proved that a 
spin-disordered phase is realized near the Mott transi- 
tion [§-fl"ll]. which could be related to the higher order 
exchange terms. As for the small A case, in addition to 
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FIG. 1: (Color online): Network of MnOe octahedrons in the 
honeycomb lattice with the definitions of interactions in the 
nearest-neighbor Ji , the second nearest-neighbor J2 , the third 
nearest-neighbor J3, and the fourth nearest neighbor J4. 



the higher order terms, the exchange pathways through 
the oxygen p state may give unexpectedly long ranged 
exchange terms and may affect the spin disordering. 

Very recently, a spin-disordered ground state is re- 
ported in a layered Mn oxide Bi3Mn40i2(N03) with 
a Mn honeycomb lattice in which the exchange inter- 
action between second neighbor Mn sites introduces a 
kind of frustration in the honeycomb lattice [r3 |. In 
this material, there is a network of MnOg octahedrons 
and Fig. [1] shows the definitions of interactions in the 
nearest-neighbor Ji, the second nearest-neighbor J2, the 
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third nearest-neighbor J3, and the fourth nearest neigh- 
bor J4. The exchange interaction J 2 between second 
neighbor sites could be derived from the Mn-O-O-Mn 
exchange pathways which becomes important in small A 
systems. The values of these interactions were studied 
both theoretically II 311 and experimentally by inelastic 
neutron scattering [14] . In these studies, only J\ and J2 
are considered by assuming that the values of J3 and J4 
are small enough to be neglected. However, since the ex- 
change interaction J 3 between the third neighbor sites is 
also derived from the Mn-O-O-Mn exchange pathways, it 
is not a trivial question whether J3 is negligible compared 
to J2 or not. In this context, it is very interesting and 
important to study the electronic structure of the Mn 
oxide, especially the values of interactions J\ - J4, us- 
ing spectroscopic methods and to reveal the origin of the 
spin-disordered state from the electronic structural view- 
point. In this paper, we investigated the electronic struc- 
ture of this material by x-ray absorption spectroscopy 
(XAS), and also estimated the values of magnetic inter- 
actions. Unrestricted Hartree-Fock calculations and a 
perturbation method revealed that the nearest neighbor 
and fourth nearest neighbor J\ and J4 are positive (anti- 
ferromagnetic), and the next nearest neighbor and third 
nearest neighbor J 2 and J3 are negative (ferromagnetic). 
In the present analysis, the ferromagnetic J2 does not 
introduce magnetic frustration in the honeycomb lattice. 
We conclude that the ferromagnetic J3 is the origin of 
magnetic frustration and the absence of long-range or- 
dering in this material. 

The synthesis of Bi3Mn40i 2 (N03) polycrystalline 
powder is described in Ref. [l2j . X-ray absorption exper- 
iments were performed at 11ID-1 (SGM) of the Cana- 
dian Light Source. The spectra were measured in the 
total-electron-yield (TEY) mode. The total energy reso- 
lution was set to 100 meV. All the spectra were measured 
at room temperature. The obtained spectrum is ana- 
lyzed by standard cluster-model calculations 15] to ob- 
tain electronic parameters. The parameters in this model 
are 3c? - 3d and 3d - 2p Coulomb interactions (Udd and 
Udc, respectively), charge-transfer energy from O 2p to 
Mn 3d states A, hopping integrals between Mn 3d and 
O 2p molecular states [V(t 2g ) and V{e g )\, and crystal 
field parameter IQDq. The superexchange interactions 
are evaluated using unrestricted Hartree-Fock calculation 
with a multi-band d — p Hamiltonian with Mn 3d and O 
2p states [l6| • The Hamiltonian is given by 



H — Hp + Hd + H p d, 



H P=Y ^kPt^Vk.la + V k P U'Ptl<yPk,l'a + H.C., 



H d - ^2 e dd^ ma d ktmc7 

k,ma 

+ J2 V k d m m-di m Jk,n-. + B.C. 



+ ( u '-f) Y d t m * d i,™* d t m 'a d h™'° 

H pd = V kj m d t ma PkM + H.c. 

k,ml,(j 

Here, d^ mcr are creation operators for the Mn 3d elec- 
trons at site i. d^ ma and p\ la are creation operators for 
Bloch electrons with momentum k which are constructed 
from the m-th component of the Mn 3c? orbitals and from 
the Z-th component of the O 2p orbitals, respectively. 
The intra-atomic Coulomb interaction between the Mn 
3d electrons is expressed using Kanamori parameters, u, 
u' , j and j' satisfying the relations u = u + j + j' and 
j = j' . The transfer integrals between the Mn 3d and O 
2p orbitals V^ m are given in terms of Slater-Koster pa- 
rameters (pda) and (pdn). The parameters determined 
by the cluster-model calculation are used as input of the 
unrestricted Hartree-Fock analysis. 

Figure [5] shows the Mn 2p XAS spectrum of 
Bi3Mn40i2(N03). There are two structures, Mn 2p 3 / 2 
— > 3d absorption at 640 - 650 eV and Mn 2pi/ 2 — > 3d 
absorption at 650 - 660 eV. The experimental spectrum 
has a sharp peak at ~ 641.5 eV characteristic of Mn 4+ , 
concluded by comparing with the reference data of Mn 2+ 
(MnO), Mn 3 + (LaMn0 3 ), and Mn 4 + (EuCoo.5Mno.5O3 
and SrMnOs) from Ref. [17] ■ This indicates that the va- 
lence of Mn is 4+ in BisMnjO^NOs), consistent with 
the valance state of Bi 3 !+ Mn4 + 0!j ! 2~(N03)~ obtained in 
Ref. [l2|. We performed configuration-interaction (CI) 
cluster-model calculations [lj| to obtain electronic pa- 
rameters. Here we fixed the following values Udd = 6.0 
eV, U dc = 7.5 eV, V(e g ) = 3.0 eV, and lODq = 1.3 eV 
and changed the value of A from 0.0 eV to 4.0 eV, as 
shown in Fig. [2] (a). The calculated spectra do not de- 
pend on the value of A so much, but the small A values of 
1.0 ± 1.0 eV reproduce the experiment most successfully 
from a peak at ~ 641.5 eV and Mn 2px/ 2 structures. 
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FIG. 2: (Color online): Mn 2p XAS spectra 
Bi3Mn40i2(N03) and comparison with the CI theory. 
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FIG. 3: (Color online): O Is XAS spectra of 
Bi3Mn40i2(NOa) together with the calculated O Is partial 
density of states by unrestricted Hartree-Fock calculations. 



for Mn 4+ specifically A = — e p + 3U. 



With this parameter set, the lowest energy state is 
found to be the conventional antiferromagnetic state 
where the first neighboring sites are antiferromagneti- 
cally coupled. In the antiferromagnetic ground state, the 
second, third, and forth neighboring sits are ferromagnet- 
ically, antiferromagnetically, and antiferromagnetically 
coupled, respectively. We have calculated energies of fer- 
romagnetic state as well as modified antiferromagnetic 
states which are obtained by exchanging spins of some 
neighboring sites in the ground state. As expected, the 
ferromagnetic state is very much higher in energy than 
the antiferromagnetic ground state. However, some of 
the modified antiferromagnetic states were found to have 
energies very close to that of the ground state. By map- 
ping the Hartree-Fock energies to the Heisenberg model 
with Ji, J2, J3, and J4, the obtained values are: 



Ji 

■h 

Ji 



9.15 meV 
-5.32 meV 
-4.80 meV 
5.77 meV 



(1) 



Figure [3] shows the O Is XAS spectrum of 
Bi3Mn40i2(N03) and the calculated O Is partial den- 
sity of states by unrestricted Hartree-Fock calculations. 
There are two structures in the O 2p - Mn 3d hybridized 
states, assigned as majority-spin e g states (e fl f) and 
minority-spin e g states (e g i), and another structure in 
the O 2p - Bi Qp states as shown in Fig. [3] There is a 
higher intensity at e 9 j_ states in the experiment than in 
the calculation, indicating that there is a substantial Bi 
6p contribution also in this structure. These assignments 
are consistent with the Mn 4+ (d 3 ) state, where majority- 
spin ti g states are occupied by electrons. 

In order to obtain the values of magnetic interactions, 
we have performed unrestricted Hartree-Fock analysis on 
the multiband d — p Hamiltonian with the Mn 3d and O 
2p orbitals. A, U, and (pda) were set to be 1.0, 6.0, and 
— 1.8 eV, respectively, on the basis of the cluster model 
analysis. The ratio ipdcr) / (pdn) is —2.16. Remaining 
transfer integrals expressed by (ppa), (ppn), (dda), and 
(ddii) are fixed at —0.6, 0.15, —0.3, and 0.15 eV, respec- 
tively, for the honeycomb lattice with the regular MnC>6 
octahedron. Here A denotes the charge-transfer energy 



While the antiferromagnetic J\ and ferromagnetic J2, 
and antiferromagnetic J4 are consistent with the antifer- 
romagnetic ground state, the ferromagnetic J3 can intro- 
duce frustration effect on it. Since the magnitude of the 
ferromagnetic J3 is comparable to those of the antiferro- 
magnetic J\ and ferromagnetic J2, the ferromagnetic J3 
by the Mn-O-O-Mn superexchange pathways is responsi- 
ble for the absence of long-range ordering in the present 
honeycomb system. 



Let us also examine the sign of J3 by considering su- 
perexchange pathways in Bi3Mii40i2(N03) using a per- 
turbation method based on the electronic structure pa- 
rameters obtained from the analysis of the Mn 2p XAS 
spectrum. Here, we will use Slater-Koster parameters, 
(ppa), (pda), (dda), and so on [l8j . 
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Jj, is given as 



J 3 



(pdir) 



[A - {ppa)/2 + (pptt)/2] 2 

[A+(ppa)/2-(ppn)/2]*) \A 
2(pdir) 4 



[A - (ppa)/2 + (pptt)/2][A + (ppa)/2 - fep»r)/2] 
11 1 



A + Up - jp 



{pda) 2 {pdTty 



[A - (ppa)/2 - (pp7r)/2P 
{pdafjpdTi) 2 

A + (ppcr)/2 + (pptt)/2] 2 
111 



u — 3j u — 2j 

2(pda) 2 (pdTr) 2 



[A - (ppa)/2 - (pp7r)/2][A + (ppa)/2 + fep»r)/2] 
11 1 1 



A 



2(pdir) 



u - 3j A + Up - j p u- 2j 
[{ppa)/2-{ppn)/2] 2 \ 1 



A 2 V A2 
4(pd7r) 4 [(ppcr)/2 - (pptt)/2] 2 



A- 



A 4 

2{pda) 2 ipdTr) 2 



Jp 



1 - 



u (A + u p ) 2 
[(ppa)/2 + (p P 7r)/2} 2 \ 1 



A 2 V A 2 J A 

A{pda) 2 {pdn) 2 [{ppa)/2 + (ppn)/2} 2 



A 2 



J 



Jp 



u 2 (A + u p y 



(2) 



Here we considered the molecular orbitals made from two 
oxygen sites. The first and second terms proportional to 
(pdir)' 1 are given by the pathway from the Mn t2 g states 
to the O 2p molecular orbitals to the Mn t^g states. The 
third and fourth terms proportional to (pder) 2 (pdir) 2 are 
given by the pathway from the Mn e g states to the O 
2p molecular orbitals to the Mn ti g states. J3 is dom- 
x inated by the negative third term and causes ferromag- 
netic interactions, which is consistent with the result of 
the unrestricted Hartree Fock calculations. 

From the neutron measurements, it was found that the 
values of interlayer interactions (J c ) are also comparable 
to Ji [3]. This is consistent with our result because J c 
is also determined by Mn-O-O-Mn pathways. 

We investigated the electronic structure of 
Bi 3 Mn 4 0i2(N03) by XAS. The valence of Mn was 
x determined to be 4+, and from CI theory we found 
that a charge-transfer energy is small in this material. 
Then we estimated the values of Ji.2,3,4 by unrestricted 
Hartree-Fock calculations and a perturbation method. 
We found antiferromagnetic J\ and ferromagnetic J2 and 
J3, leading to the existence of magnetic frustration and 
the absence of long-range ordering, as experimentally 
confirmed. 
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